arXiv:1503.07592v2 [hep-ph] 7 Jul 2015 


TIFR/TH/15-08, LPSC-15/081 


Probing the NMSSM via Higgs boson signatures from stop 

cascade decays at the LHC 


Amit Chakraborty*^, Dilip Kumar Ghosh*’, Subhadeep MondaP, 

Sujoy Poddar'*, Dipan Sengupta® 

® Department of Theoretical Physics, Tata Institute of Fundamental Research, 

1, Horn! Bhabha Road, Mumbai-400005, India. 

*’ Department of Theoretical Physics, Indian Association for the Cultivation of Science, 

2A & 2B, Raja S.C.Mullick Road, Jadavpur, Kolkata 700 032, India. 

® Harish-Chandra Research Institute, Chhatnag Road, Jhusi, Allahabad 211019, India. 

'*Netaji Nagar Day College, 170/436 N.S.C. Bose Road, Kolkata 700092, India. 

® Laboratoire de Physique Subatomique et de Cosmologie, Universite Grenoble-Alpes, CNRS/IN2P3, 
53 Avenue des Martyrs, F-38026 Grenoble Cedex, France. 


Abstract 

Higgs signatures from the cascade decays of light stops are an interesting possibility in the next to minimal 
supersymmetric standard model (NMSSM). We investigate the potential reach of the light stop mass at the 
13 TeV run of the LHC by means of five NMSSM benchmark points where this signature is dominant. These 
benchmark points are compatible with current Higgs coupling measurements, LHC constraints, dark matter 
relic density and direct detection constraints. We consider single and di-lepton search strategies, as well as 
the jet-substructure technique to reconstruct the Higgs bosons. We find that one can probe stop masses 
up to 1.2 TeV with 300 fb“* luminosity via the di-lepton channel, while with the jet-substructure method, 
stop masses up to 1 TeV can be probed with 300 fb~* luminosity. We also investigate the possibility of the 
appearance of multiple Higgs peaks over the background in the fat-jet mass distribution, and conclude that 
such a possibility is viable only at the high luminosity run of 13 TeV LHC. 


1 Introduction 

The discovery of a Higgs-like particle with mass close to 125 GeV by both the ATLAS [1] and CMS [2] 
collaborations has led us to a new crossroad. Measurements of its spin, CP quantum number and couplings 
to the standard model (SM) particles so far are consistent with SM predictions for the Higgs boson [3,4]. 
The natural question is whether this Higgs-like particle is also one of the CP-even Higgs bosons of the 
supersymmetric (SUSY) models (for a review of supersymmetry, see Refs. [5-7]). The primary goal of the 
run-H of the LHC therefore is to look for new particles as well as investigate the properties and couplings of 
this Higgs boson more precisely. 

In the minimal supersymmetric standard model (MSSM), in order to push the light CP-even Higgs boson 
mass up to 125 GeV, a sizeable amount of quantum correction is required. The dominant contribution comes 
from the light third generation squarks, namely the stops. However, to achieve a 125 GeV Higgs boson in the 
MSSM one requires to have maximally mixed or heavy stops (close to 1 TeV). This requirement, along with 
the lower limits imposed on the different SUSY particles (specially the gluino and the light stop) after the 8 
TeV run at the LHC have pushed the fine-tuning problem of the SM to an uncomfortable level [8] . The next 
to minimal supersymmetric standard model (NMSSM) [9] provides a solution to this problem. The model 
extends the MSSM field content by the addition of a gauge singlet chiral superfield (S). As a consequence 
of the coupling of this new gauge singlet superfield S with the Higgs doublet superfields Hu and Hd, the 
AHuHdS term in the NMSSM superpotential, the tree level Higgs boson mass receives an extra tree level 


contribution which is proportional to the square of the singlet-doublet coupling A. Therefore a relatively 
smaller radiative correction is required in the NMSSM to obtain a 125 GeV Higgs. This ensures that the 
fine-tuning issue is diluted to some extent in the NMSSM compared to the MSSM [10-15]. 

The phenomenology of the NMSSM is much richer than MSSM due to the presence of the additional 
superfield S. The Higgs sector of the NMSSM consists of five neutral (three CP-even and two CP-odd 

H- 1 . 2 ) and two charged scalars. The neutral component of the scalar part of the gauge singlet superfield 

S mixes with the neutral components of the Higgs fields, and thus there can be a significant doublet-singlet 
mixing after the Higgs mass matrix is diagonalized. For certain choices of the model parameters, one can 
find a light Higgs boson which is predominantly singlet like. Therefore regardless of the CP properties, its 
coupling with the SM gauge bosons will be highly suppressed. Therefore the constraints imposed by LEP 
cannot exclude such a singlet like light Higgs boson. Thus having additional Higgses in the vicinity of the 
SM-like 125 GeV Higgs boson is an attractive possibility in the NMSSM [16-21]. While the spin, parity and 
coupling measurements of the discovered Higgs boson indicate that it is the SM-like Higgs, it could very well 
be the second lightest CP-even Higgs boson of the NMSSM, with SM like couplings to fermions and gauge 
bosons. 

The Higgs bosons in the NMSSM can arise from the cascade decays of stops and heavy neutralinos. A 
wide range of papers have studied the production of Higgs bosons and neutralino cascade decays to Higgs 
bosons in the NMSSM, and we refer some of them in [22-28]. The prospect of having additional NMSSM 
Higgs bosons in the vicinity of the 125 GeV Higgs boson was considered in [11], where the authors carefully 
analyzed the interplay of NMSSM Higgs bosons and the stop sector and concluded that one could achieve 
a modest level of fine tuning with a light stop below 1 TeV. Moreover, the possibility of the production of 
multiple Higgses in NMSSM was also considered in [29], where the authors considered pair production of 
SM like Higgs bosons from squark and gluino cascade decays. 

In a R-parity conserved SUSY scenario, the lightest neutral, stable SUSY particle (LSP) serves as a good 
cold dark matter (DM) candidate. The presence of a singlet superfield in the NMSSM, that gives rise to an 
extra singlino component in the neutralino sector, results in a very interesting DM phenomenology. Moreover 
unlike the MSSM, a very light DM matter candidate as favored by some of the DM experiments [30,31] can 
also be accommodated in NMSSM. Therefore a comprehensive phenomenological study of the interplay of 
the Higgs sector and the neutralino sector in the NMSSM is of prime importance at the LHC. 

At this point it is pertinent to discuss some of the LHC analyses dedicated to stop searches at the 8 
TeV run of LHC relevant to this study. Most of these searches by the ATLAS and CMS collaborations 
have been conducted in the framework of simplified topologies. The CMS collaboration has studied the 
prospect of discovering the light ti assuming 100% branching ratio (BR) for ti —>■ ty?, and placed the limit 
on the stop mass up to 650 GeV for massless neutralinos in the single lepton -|- jets + channel [32]. The 
decay of the heavier stop (t 2 ) to ti and the SM-like Higgs (t 2 —>■ tiHsM) has also been studied by the CMS 
collaboration [33]. The limits at 8 TeV are rather weak, with the heavier stop up to 575 GeV and the lighter 
stop up to 400 GeV being ruled out with 19.5 fb“^ of integrated luminosity. Additionally, in the framework of 
gauge mediated supersymmetry breaking (GMSB), the channel ti —>■ ty?, followed by the decay of Xi to Hsm 
and a Gravitino (G) was studied in [34]. The limits placed on the light stop mass for this particular scenario 
is about 400 GeV. Furthermore, the pair production of chargino and neutralino with the decay X 2 XiHsm 
has also been studied by the ATLAS collaboration, assuming degenerate masses for xt X 2 - The limit 
on m^o was set to 250 GeV [35]. It has to be remembered that most of these experimental studies have 
been performed assuming that the branching ratios for all the topologies are 100%. A few phenomenological 
studies with boosted top and Higgs have also been performed for 14 TeV run of LHC in the simplified model 
scenarios and in the context of some specific models as well. In particular, the prospects of x^X 2 HsmWx? 
with boosted Higgs was studied in [36] in the context of the constrained minimal supersymmetric standard 
model (CMSSM). It was concluded that xtiX 2 masses up to 400 GeV could be probed at 14 TeV LHC in 
this framework. The decay of heavier stops to t 2 —>■ UHsm was considered in [37]. The author concluded 
that (2 masses up to 1 TeV could be probed at 14 TeV with 100 fb“^ luminosity, in a simplified model 
scenario. 

In this paper, we consider the production of Higgs bosons from the cascade decay of the ti in the NMSSM 
framework. We investigate the prospects of light stop pair production followed by the decay of ti to a top 
quark and any one of the heavier neutralinos {Xi =2 s): which subsequently decays to one of the neutral 
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Higgs states (Hi,Hsm = H 2 , Ai)^ and the lightest neutralino, i.e., x? Xi + Hi_sm/Ai. The produced 
Higgs in this case will be boosted if there is a significant mass gap between m^o and m^o. Here we consider 
five benchmark points (BP) compatible with the current Higgs coupling measurements, DM relic density 
constraint, and direct detection constraints. The main characteristic feature of all these BPs is that the 
decay of the heavier neutralinos to the Higgses is dominant. We apply the jet substructure technique to 
reconstruct these boosted Higgs bosons to estimate the reach of stop searches in the context of 13 TeV run 
of LHC. We also analyze the possibility of the appearance of multiple Higgs peaks over the background in 
the fat-jet mass distribution. It is observed that the appearance of such peaks are only viable at very high 
luminosity run of LHC. However it is worth pursuing such a possibility to reveal the presence physics beyond 
the SM. We also compare our results with the generic lepton(s), jets + missing energy searches at the LHC. 

The rest of the paper is organised as follows. In section 2 we discuss the choice of benchmark points, the 
relevant constraints imposed and briefly review some of the literature relevant to this study. In section 3, we 
discuss our collider analysis, while in section 4 we present our results and finally we conclude in section 5. 


2 Benchmark points and constraints 

As already discussed, we focus on the production of the Higgs bosons from the cascade decay of the ti. Thus 
we start with ti pair production followed by the cascade decay of the stop to neutralinos and Higgs. The 
final state therefore contains top quarks, Higgs bosons, and LSPs. The Higgs decays dominantly to pair of 
b-quarks, while we let the top quark decay inclusively. We present a sample Feynman diagram of the process 
of our interest in Fig. I. Few representative cascade decay modes are shown below, where X collectively 
denotes the decay products of the top quarks and also other hnal state particles except the Higgses: 
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Figure I: A sample Feynman diagram for the channel under consideration. 


The decay of the next-to lightest supersymmetric particle (NLSP) to LSP -b X (X = Z/Hsm) in the 
NMSSM framework has been already addressed in [38] , where the authors pointed out that for a singlino like 
LSP satisfying DM constraints, this decay would lead to a reduction of missing transverse energy making it 
difficult to probe at the LHC. The authors also looked at the prospects of Higgs production from the cascade 
decay of squarks and noted that small values of A and k were required to satisfy the invisible Higgs decay 

^Throughout this paper we will denote the SM-like Higgs as HgM = H 2 and the rest of the Higgs spectrum contains 

H=*=, Hi, 3 , Ai, 2 . 
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constraint by suppressing the mixing of the singlet like Xi with Higgsinos, for cases where m^o ~ 5 GeV. 
The authors showed that the jjbb background was significantly large at 13 TeV LHC and so they proposed 
that tagging multiple b-jets along with the requirement of a large number of jets can be useful to reduce 
the large background. The decay X 2,3 XiHsm was also discussed in [39], where the authors also argued 
along the same lines. The authors looked at pp —>■ with subsequent decay of the neutalino states 

to the SM-like Higgs. They showed that four b-tagged jets in addition to can be a viable channel at 14 
TeV LHC. The authors also suggested ways of distinguishing MSSM from NMSSM by observing dark matter 
patterns and neutralino branching ratios at the colliders. This particular region of parameter space specific 
to the NMSSM has therefore received significant attention. From the above studies it is also evident that 
probing such a channel at the LHC environment is a very difficult task. A detailed collider study at 13 TeV 
LHC is preformed here. 

The benchmark points for this study are chosen with the intention of probing the regions of parameter 
space where the decay of X 2 3 LSP + Higgs is dominant. We present our five representative benchmark 
points in Table 1. The branching ratios of the relevant decay modes and some of the crucial observables 
corresponding to the benchmark points are tabulated in Table 2. The particle spectrum and decays are 
generated using NMSSMTOOLS-4.2.1 [40]. Since we are interested in the decay of the light stop, we decouple 
the rest of the spectrum, i.e, we set the mass of the gluino and the rest of the squarks and sleptons at 3 TeV. 
Furthermore, we want the light stop to be predominantly left handed and therefore we set the soft masses 
for the right handed 3rd generation squarks to be at 3 TeV, so that the decays ti —tX 2,3 ti —hxi 
are dominant. It is pertinent to note that a detailed collider analysis for the left handed light stop was 
performed by some of the same authors in an earlier work using the boosted top technique in the context of 
phenomenological MSSM [41]. Since our primary channel of interest for the subsequent stage of the cascade 
is X 2,3 XiHsm/Hi/Ai, we allow top to decay inclusively, while Higgs decays dominantly to bb. It has to 
be remembered that the whole scenario is additionally complicated by the fact that one has to satisfy the 
DM relic density and also DM direct detection constraints, if Xi is assumed to be a valid DM candidate. 

In order to ensure that X 23 ii®'S a large higgsino component, we fix Mi and M 2 to 1.5 TeV for all 
benchmark points and vary /ieS in the range of 260-280 GeV. This also ensures that X 4 5 ^re predominantly 
gaugino like and have mass close to 1.5 TeV. Thus they play no role in the cascade decay of the light stop. 
On the other hand, the large value of A ensures that the lightest neutralino is predominantly singlino for all 
of the parameter space points. The two decay modes that compete here are X 2 3 XiZ and X 2,3 XiHi.sm, 
when the mass difference, m^o ^ — nijj,o > mHj sm is satisfied. In the MSSM, this competition is predominantly 
determined by the gaugino-higgsino components of the neutralino, such that the decay of heavier neutralino 
states to lighter ones and Higgs has full gauge strength if one of the neutralino is gaugino like while the other 
is higgsino like. However, in the NMSSM, the LSP can have a significant amount of singlino admixture, and 
this is the case for our all the benchmark points. It is known that the Z boson only couples to the higgsino 
like component of the lightest neutralino and therefore the XiX^^ coupling is suppressed for a dominantly 
singlino like neutralino. Hence X 2 XiHsm dominates in all our benchmark points (see Table 2). Here, we 
perform a small scan over the NMSSM parameter space and consider the parameters presented in Table 1 
such that the relevant decay modes are dominant and all the present low energy constraints are satisfied. In 
the Higgs sector, we ensure that the LEP constraints as well as SM Higgs coupling constraints at the end 
of 8 TeV run of the LHC are satisfied [3,4,42]. The lightest Higgs boson being dominantly singlet like, the 
HiZZ coupling is suppressed ensuring that the LEP bounds are easily satisfied. A consequence of the singlet 
like Hi is that the branching ratio of X 2 3 XiHi is rather small (see Table 2). We also note that Ma, 
the pseudoscalar mass parameter, plays an important role to satisfy the DM constraints. For the parameter 
space of our interest, the LSP annihilates dominantly via the pseudoscalar, and therefore a small tuning in 
the mass of Ma is required to obtain the DM relic density within 2a of WMAP/PLANCK data [43,44]. 
Additionally, we make sure that the benchmark points of our choice satisfies the direct detection cross-section 
bounds obtained from the LUX [45] collaboration. 

All the benchmark points satisfy the following set of constraints: 

• The mass of the SM-like Higgs must be within the range 123 GeV < Mhsm < 129 GeV. Due to the 
small difference in the central values of the ATLAS and CMS measurements and also considering the 
theoretical uncertainty in Higgs mass calculation, we use 126 ± 3 GeV as a conservative estimate [46,47]. 

• The SM-like Higgs couplings to be within 2a of the measured values provided by the ATLAS and CMS 
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PI 

P2 

P3 

P4 

P5 

™Q3 

830 

920 

1000 

1180 

1350 

tan P 

10.2 

12.2 

10.4 

10.3 

10.2 

A 

0.34 

0.15 

0.52 

0.6 

0.44 

K, 

0.06 

0.03 

0.09 

0.09 

0.07 

Aa 

2700 

2700 

2700 

2700 

2700 

A« 

-1200 

-1200 

-1200 

-1200 

-1200 

Pe// 

280 

264 

262 

268 

277 


124.3 

124.9 

124.3 

125.1 

124 

rriHi 

67.1 

67.1 

62.3 

73.8 

67.2 

ruA-i 

214 

226.8 

194.4 

154.0 

193 

mr 

804.2 

908.2 

1003.7 

1211 

1392.5 

IHk 

bi 

821.4 

923.1 

1018 

1226 

1408.6 

Xi 

98.6 

107.2 

87.3 

77.6 

87.4 


290.1 

268.0 

282.7 

293.2 

292.4 

m^o 

293.5 

273 

283.1 

294.2 

295.1 

m ± 

Xi 

284.7 

268.5 

266.5 

273.2 

282.5 


Table 1: The parameters and masses for the five benchmark points. All the other parameters are set to their 
fixed values as described in the text. All masses are in GeV unit. 



PI 

P2 

P3 

P4 

P5 


0.95, 0.96 

0.94, 0.94 

0.96, 0.96 

0.92, 0.91 

0.92, 0.91 

glliggF), g.ll{VBFIVH) 

1.0, 0.99 

0.99, 0.97 

0.99, 0.98 

1.0, 1.0 

1.0, 1.0 

giliUH), ^1\{VBF/VH) 

1.0, 0.97 

0.95, 0.95 

1.0, 1.0 

1.01, 0.93 

1.0, 1.0 

Q,hP, asi{xW'^^, cm^) 

0.123, 6.4 

0.122, 1.9 

0.115, 8.1 

0.125, 1.8 

0.126, 9.3 

BR(ti ^ty§,3,4)(%) 

92 

94 

97 

81 

88 

BR(ti ^bx^2)(%) 

7 

6 

3 

4 

4 

BR(x^^X?Hi,sm)(%) 

6, 60 

8, 61 

20, 57 

18, 49 

14, 54 

BR(x§^X?Hi,sm)(%) 

3, 14 

8, 1 

3, 15 

6, 16 

4, 16.3 

BR(Hi ^ bb)(%) 

86 

69 

91 

90 

85 

BR(H2 ^ bb)(%) 

64 

64 

64 

64 

66 

BR(Ai ^ bb)(%) 

1 

1 

1 

90 

1.5 


Table 2: Some of the Higgs signal strengths, dark matter relic density, direct detection constraints and the 
relevant branching ratios for the particular study. 


collaborations at the end of 7 + 8 TeV run of LHC with approx. 25 fb ^ of data [3,4]( See Table 3 for 
more details.) 

• The HR of the rare b-decays B —>■ Xgy and Bg —?► within 2a of the experimental results [53]. 

We assume: 
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Signal Strength (/r) 

Channel 

CMS 

ATLAS 

77 

1.12 ±0.24 [4] 

1.17 ±0.27 [48] 

WW 

0.83 ±0.21 [4] 

1.09l°J? [49] 

ZZ 

1.0 ± 0.29 [4] 

1.44l° j° [50] 

bb 

0.84 ± 0.44 [4] 

0.52 ±0.32 ±0.24 [51] 

TT 

0.91 ±0.28 [4] 

1.43t°;« [52] 


Table 3: Updated results on Higgs coupling measurements by the ATLAS and CMS collaborations at the 
end of 7 + 8 TeV run of the LHC with approx. 25 fb“^ of data. 


2.77 X 10-'‘ < BR(B ^ X^y) < 4.09 x 10"^ 

1.0 X 10"® < BR(Bs Ai+Ai”) < 5.2 x 10"® 

• The LEP lower bound on the chargino mass, M_^± > 103.5 GeV [42]. 

• We apply the upper bounds on the DM direct detection cross-section from LUX [45] and allow the DM 
relic density to be within 2a of the WMAP/PLANCK [43,44] data. Note that, precise value of DM 
relic density do not have any significant impact on our final results. We assume. 


0.115 < HDMh^ < 0.126. 


Before we end this section, we would like to remind the readers that for all our benchmark points, the 
branching ratios for the decay ti —>■ tx 2 3 ^.re significantly large owing to the large higgsino component in 
X 2,3 (see Table 2). The dominantly Higgsino-like X 2 ^^lel largely singlino-like Xi ensure that X 2 XiHsm is 
significantly large, which is necessary for our collider strategy. Since Hsm is SM-like, it has a large branching 
ratio to bb, while Hi, due to the large singlet component and reduced coupling to ZZ also dominantly decays 
to bb. From Table 2, it is evident that all the reduced signal strengths of the SM-like Higgs are well within 
2tT of the limits imposed from the Higgs coupling measurements at the LHC. With these comments on the 
benchmark points, we now proceed to discuss our signal and background collider analysis in the next section. 


3 Analysis of the signal and background 

In this section we describe the strategy for the collider study. The entire analysis is divided into two parts: 
first we use the jet substructure technique to analyze the boosted Higgs scenario. This is followed by the 
standard leptonic search strategies at the LHC. We then compare these two methods and discuss the merits 
and demerits of each of these two procedures. 

3.1 Jet substructure analysis 

The Higgs tagging technique, introduced in the work of Butterworth-Davison-Rubin-Salam (BDRS) [54], 
has emerged as a powerful tool to study physics in the boosted regime. The efficacy of the BDRS tagger 
lies in the sharp resolution of the Higgs peak in the H —>■ bb channel. In the analysis, we apply the BDRS 
Higgs tagger to tag the Higgs originating from the decay of the heavier neutralinos. From Table 1, it is clear 
that the mass difference m^o — m^o is not significantly big enough to generate a large boost to the Higgs. 
However, being part of a cascade decay, the boost to the Higgs boson originate from the two step decay 
chain of the stop. To demonstrate the larger boost in the signal process as compared to the background. 
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Figure 2: px-distribution of the leading two fat-jets for 13 TeV the LHC. We reconstruct the jets using the 
C/A algorithm with jet radius R = 1.2. The benchmark points BPl and BP5 are as tabulated in Table 2. 


we present, in Fig. 2, the px distribution of the leading two “fat-jets” reconstructed using the Cambridge- 
Aachen (C/A) [55] jet algorithm with the jet radius R=1.2, for two benchmark points BPl and BP5, as well 
as the SM tf -I- jets background. One can observe that for the signal, the peak of the distribution is over 400 
GeV, while the background falls sharply from 200 GeV. We notice that with the increase in the stop mass 
(BPl to BP5), the fat-jet distribution moves to the higher px region, as expected. To optimize the Higgs 
tagging efficiency we considered the values of the fat-jet radius R to be 1, 1.2 and 1.5 and concluded that 
R=1.2 is the best choice in terms of signal to background ratio. Furthermore we varied the jet px threshold, 
considering three different values 200 GeV, 400 GeV, 600 GeV and arrived at the conclusion that a jet radius 
R=1.2 and a jet px threshold of 200 GeV gives the best efficiency in terms of enhancement in the signal 
and reduction of background. Since the boost to the Higgs boson in our case is not significantly large, we 
achieve a modest Higgs tagging efficiency (e ^ 5-6 % for px ~ 200 - 300 GeV). However this is an order of 
magnitude larger compared to the largest background in terms of cross section, namely tt-fjets (see Table 4 
and 5). An added advantage of using the BDRS Higgs tagger is that it uses a filtering technique to discard 
the contamination from the underlying events, which in turn helps to improve the resolution of the Higgs 
peak. As a consequence of the superior mass resolution, it is also a powerful tool against the backgrounds 
containing multiple b-jets like tt -I- jets, Z(—bb) -I- jets, ttbb, ttW, ttZ. The SM backgrounds that we 
consider are listed below, 

• tt -I- jets : The decay of top quark yield two b-jets in the final state and thus the reconstructed fat-jet 
could in principle fake a Higgs signal. 

• Z(— 7 bb) -I- jets: In this case although the fat-jet mass distribution should peak at the Z boson mass 
there might be a large tail. Owing to the large Z -|- jets cross section, this could potentially be significant 
in the end. 

• ttbb, ttZ, ttW : These backgrounds contain at least 2 b-jets, which may be tagged as a fake SM Higgs 
boson. 

• ttHsM : This constitutes the irreducible background to our signal with at least one tagged Higgs jet 
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and ^rp. However the Higgs in this case is not expected to be as boosted as the signal. 

We use NMSSMTOOLS-4.2.1 [40,56-58] to generate the NMSSM mass spectrum, while signal and back¬ 
ground events are simulated using MADGRAPH5 [59] and passed to PYTHIA6 (version 6.4.26) [60] for 
showering and hadronization. For the signal (titj -I- 2 jets) and tt/Z/W -I- jets (up to 2 jets), we perform the 
matrix element parton shower (ME-PS) merging using the MLM prescription [61] with a merging parameter 
of 30 GeV. We use Fast Jet-3.1.0 [62] for jet reconstruction, and implementation of our jet substructure 
analysis. The rare b-decay observables, DM relic density and direct detection cross-sections are calculated 
using NMSSMTOOLS. 

We observe that even after tagging of the Higgs jet, there is still a significant amount of hadronic activity 
in the signal process. This activity can be used as an effective handle against the background. Therefore 
the jet reconstruction procedure in our case is a two step process. We first reconstruct fat-jets using the 
Gambridge-Aachen (C/A) [55] algorithm using a jet radius of R=1.2 with a jet px threshold of 200 GeV 
followed by the use of the BDRS [54] Higgs tagger with a default mass drop criteria of ^ = 0.67 and the 
filtering parameter Rgit = min(R|35/2, 0.3). Once one (or multiple) SM-like Higgs boson(s) are tagged, we 
remove those particles from the fat-jet list and recluster the remaining stable hadrons using anti-kx [63] jet 
algorithm with a jet radius of R=0.4, px > 20 GeV. We select jets with plp > 30 GeV and pseudorapidity 
|?7j| <3 for further collider analysis. Leptons (electrons and muons) are selected with pi^ > 10 GeV and 
|? 7 ^| < 2.5. Isolation of leptons are performed by demanding that the sum of the scalar px of all stable visible 
particles within a cone of radius AR = 0.2 around the lepton should not exceed 10% (15% for muons) of 
Px (Px)- missing transverse energy is constructed using all the stable final state visible particles with 
all jets with plp > 20 GeV and |? 7 j| < 4.5, and all leptons with p^ > 10 GeV and |? 7 ^| < 2.5. Following the 
ATLAS collaborations [64], we consider a px-dependent b-tagging efficiency : 


0 

for 

p^J. < 20 GeV 

0.6 

for 

20 GeV < p!(. < 50 GeV 

0.75 

for 

50 GeV < p!^ < 400 GeV 

0.5 

for 

PX > 400 GeV 


After generating events, the following set of kinematic cuts are applied, 

• Cl: We demand that at least one of the fat-jets is identified as the SM-like Higgs boson with the criteria 
110 GeV < mj < 140 GeV, where mj denotes the fat-jet mass. The fat-jet is deemed to be a Higgs-like 
jet if and only if two b-jets are identified inside the fat-jet. In the signal process, however, there are 
additional Higgs bosons which can be tagged and therefore at the end of all selection cuts we look 
at the mj distribution, to find if additional peaks are visible over the background. The Higgs boson 
tagging ensures that most of the backgrounds with multiple b-jets in the final state are significantly 
reduced. 

• C2: In this work we consider a hadronic final state, and therefore we impose a lepton veto. 

• C3: For the signal, even after at least one SM-like Higgs boson is tagged, there are significant number 
of b-jets originating from the top quarks. However, these extra b-jets are unlikely to be present in the 
backgrounds like Z (—>• bb) -I- jets, tt -I- jets after satisfying the Higgs mass criteria. Therefore, we 
demand at least one additional b-jet in the event. 

• C4: An additional selection criteria of Hx = Sjplp > 800 GeV is imposed. The jets in this case 
are the anti — kx jets constructed by reclustering the rest of the hadronic sample after removing the 
hadrons contributing to the Higgs tagged jet as mentioned earlier. The larger hadronic activity in 
the signal implies a larger value of Hx compared to the background. This cut is particularly useful in 
suppressing the irreducible ttHsM background. In the left panel of Fig. 3, we present the Hx distribution 
corresponding to two signal benchmark points BPl and BP5, along with the two dominant backgrounds 
tt -I- jets and ttHsM- 

• G5: As a final selection criteria we impose a ;^rp cut of 175 GeV. The signal jirp distribution is harder 
compared to that of the SM backgrounds because of the heavy LSP in the final state of the signal 



events. In the right panel of Fig. 3, the pirp distribution for the signal and tt + jets are presented to 
illustrate this feature. 




Figure 3: Left panel shows the distribution of Ht where Ht = Sjplp with the anti-kx jets selected after 
imposing the basic Higgs tagging criteria, while right panel displays the distribution of the pirp at the 13 TeV 


LHC. 


3.2 Alternative analysis: Final state with leptons 

In order to understand whether the jet substructure technique is more effective than other search strategies, 
we perform a separate collider analysis with leptonic final states. Here we probe two such possible scenarios, 
one dedicated to the single lepton hnal state (denoted as Case-I), while other being the di-leptonic channel 
(denoted as Case-H). The SM backgrounds for both these channels are tt +jets, W(—>■ Iv) +jets, Z +jets, 
tfZ(—>■ ttW (where the leptons can originate either from W or from the top) and the irreducible 

ttHsM- To analyze these two channels, we devise two sets of search strategies and optimize our selection cuts 
to obtain a good signal significance at the LHC. We use FASTJET to reconstruct the jets with anti-kx [63] 
jet algorithm with a jet radius of R=0.4 and p™'" > 20 GeV. We select jets with pip > 30 GeV and |?7j| <3 
for further analysis. Below, we list the optimized selection cuts for both the cases separately. 

1. Case I (Single lepton) : 

• C6: We demand at least 4 jets with pi^ > 200 GeV, pi^ > 150 GeV, piji > 100 GeV, piji > 50 GeV, 
since the signal has a significantly larger number of hard jets as compared to the background. 

• C7: Additionally, we demand at least 2 b-tagged jets, to suppress backgrounds like W (—>■ li^)+jets. 

• C8: A single isolated lepton is required with p^ > 20 GeV and \r]e\ < 2.5. 

• G9 : Since the leptonic mode is under investigation, a large ^rp is required to suppress backgrounds 
like tt + jets, ttHsM- Therefore a piq. > 400 GeV is imposed. 

• CIO : Due to the signihcantly larger jet activity in the signal, we impose Hx = Splp > 600 GeV, 
where jets satisfying the basic selection criteria are considered here. 
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Figure 4: The Mt 2 distribution at 13 TeV LHC energy for the benchmark points PI, P5 and the SM 
backgrounds tfHsM and tt + jets. 


• Cll : The SM process ttHsM, with one of the tops decaying leptonically is the largest background. 
To suppress this, we demand that the transverse mass of the lepton and missing transverse energy 

of the system, Mt(^, x x (1 — cos0(.f, ;^rp)) > 200 GeV. This ensures that the 

single lepton final states originating from the W boson are suppressed as the end point of the 
Mt(^, distribution is expected to be bounded from above by the mass of the W boson. How¬ 
ever, for the signal the presence of the additional Xi ensures that the end point of the Mt(-^,p1x) 
distribution is shifted beyond the W boson mass. 

• C12: Events with number of jets greater than six are accepted. 

2. Case II (Di-lepton) The di-lepton channel, although being cleaner compared to the single lepton 
channel, has the disadvantage of a lower signal cross section. We optimize the di-lepton channel by 
imposing the following set of cuts. 

• C13 : We demand at least 3 jets with pijl > 150 GeV, p^ > 100 GeV, p^ > 50 GeV. 

• C14 : At least two b-tagged jets are required. 

• C15 : We demand two isolated leptons with p^ > 20 GeV and \r]^\ < 2.5. 

• C16 : In order to suppress the backgrounds like tt -f jets, ttHsM, which yield di-leptons from the 
W decay, we use the variable Mx 2 [65,66], defined as, 


Mx2(l^X^l^T^yT) 


= min 


max{Mx(lt", fx), Mx(lt^x^ ^t)} 


( 1 ) 


where il and £2 are the two isolated leptons, ^x the total missing transverse momentum of 
the event and Mx(l^G 1^^) is the transverse mass of the system, defined as 
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(j) being the (azimuthal) angle between and . In the definition of Eq.l, and are the 
two hypothetical splits of the total missing transverse momentum. It is assumed that the mass 
of the invisible particles are zero [67]. The use of di-leptonic Mt 2 to suppress SM backgrounds 
has been considered previously in a host of other works, for example see Refs. [41,68-70]. For 
leptons originating from the decay of the top quark, the value of Mt 2 is bounded by the W mass, 
whereas for signal the presence of Xi ensures that the distribution is shifted beyond the W mass. 
The signal and background distribution illustrating this feature is presented in Fig. 4. In this 
work, we set Mt 2 > 150 GeV. 

• C17 : Additionally we impose a criteria on the total missing transverse momentum, jirp > 400 GeV. 


4 Results 

In this section, we summarize our findings based on the collider simulation described in the previous section. 
We begin our discussion with the the jet substructure analysis. The event summary for the signal after 
individual selection cuts are presented in Table 4, while in Table 5 the event summary for the backgrounds 
are tabulated. The second column of Table 4 corresponds to the NLO pair production cross-section of the 
light stop, as obtained from the official LHC supersymmetry cross section working group [71], corresponding 
to the benchmark points in Table 2. The cross sections for tt -I- jets are obtained using MadgraphS with 
an appropriate K-factor for 13 TeV LHC, while for the cross section of ttHsM process we use LHC Higgs 
cross-section working group report [72]. Although all the other backgrounds like W/Z-l-jets, ttbb, ttZ, ttW 
as mentioned in the previous section were simulated, we present only tt -I- jets and ttHsM, as these are the 
dominant ones. It was checked that the choices of the selection cuts are such that other processes do not 
contribute to the final background cross section. 




Effective cross-section after the cuts (in fb) 

Signal 

Production 

c.s. (fb) 

Cl 

C2 

C3 

C4 

C5 

PI 

27.2 

1.18 

0.58 

0.48 

0.25 

0.146 

P2 

12 

0.42 

0.21 

0.17 

0.11 

0.0794 

P3 

6 

0.22 

0.11 

0.09 

0.07 

0.0529 

P4 

1.5 

0.0497 

0.0262 

0.0225 

0.0193 

0.0161 

P5 

0.5 

0.0258 

0.0139 

0.0121 

0.0115 

0.0096 


Table 4: Jet Substructure: Event summary for the signal after individual cuts as described in the text. 
The simulation is performed at 13 TeV LHC energy. From 3rd column, we show the effective cross-section 
(c.s.) after all the event selection cuts have been applied. 


From the 3rd column to the 7th column (C1-C5), we present the cut flow, normalized to the cross sections 
after imposing the selection cuts, as discussed in the previous section. Similarly, in Table 5, we discuss the 
event summary for the backgrounds after each cut. As can be observed from the 3rd column of Table 
4, about 3-6% percent of the events have at least one tagged Higgs for the signal events. The tt -|- jets 
background is reduced by 99.89% (see Table 5, 3rd column) while backgrounds like W/Z-|-jets, ttbb, ttZ, 
ttW (not shown in the Table ) are reduced to negligible amounts. The largest background, as expected is 
ttHsM, with a tagging efficiency of about 2.8%. Therefore to reduce this background the subsequent cuts 
C2-C5 play an important role. From Table 5, we notice that the principal cut to suppress the rest of the tt 
-bjets background is the demand of extra b-jets (cut C3) after satisfying the criteria of at least one SM-like 
Higgs. As mentioned earlier, these extra b-jets are unlikely to be present in the background and we observe 
that about 98% of the remaining tt -bjets background is further rejected. In signal and ttHsM however these 
extra b-jets are naturally present, and therefore they are not significantly suppressed. To suppress the ttHsM 
background, we employ the cut on Ht and pirp as noted in the previous section. Even after the Higgs tag the 
presence of significantly high energetic jets from the stop decay as compared to ttHsM and tt -bjets ensures 
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Effective cross-section after the cuts (in fb) 

SM bkgs 

Production 

c.s. (fb) 

Cl 

C2 

C3 

C4 

C5 

tf-b jets 

ttHsM 

700000 

500 

751.02 

14.45 

554.85 

9.76 

7.5 

7.47 

0 

0.5 

0 

0.028 

Total 

bkg 


0.028 


Table 5: Jet Substructure; Event summary for the backgrounds (bkgs) after individual cuts as described 
in the text. Numbers in the last column show the final cross-section (c.s.) after all the event selection cuts 
are applied. For the tt-b jets background, we generate a matched sample of tf-b 0 jet, tt-b 1 jet and tt-b 2 jets 
using Madgraph. We also consider IT-bjets and Z-bjets backgrounds, however both of them are identically 
zero at the end of C5. 



Signal(Ns) ( Background(NB)) 

Significance).?) for 

K = 10% 


m*! (GeV) 

100 fb“i 

300 fb“i 

1000 fb"i 

100 fb”i 

300 fb“i 

1000 fb"i 

PI 

804.3 

14.6 (2.8) 

43.8 (8.4) 

146 (28) 

8.6 

14.5 

24.37 

P2 

908.2 

7.94 (2.8) 

23.82 (8.4) 

79.4 (28) 

2.75 

7.9 

13.25 

P3 

1003.7 

5.29 (2.8) 

15.87 (8.4) 

52.9 (28) 

1.83 

5.25 

8.83 

P4 

1211.5 

1.61 (2.8) 

4.83 (8.4) 

16.1 (28) 

0.55 

1.6 

2.69 

P5 

1392.5 

0.96 (2.8) 

2.88 (8.4) 

9.6 (28) 

0.33 

0.95 

1.6 


Table 6: Jet Substructure: The summary of our signal and backgrounds. Columns 3-5 show the number 
of signal (total background) events for three values of the integrated luminosity: 100 fb“^, 300 fb“^ and 
1000 fb“^ at 13 TeV LHC. The columns 6-8 show the statistical significance of our signal for the above three 
integrated luminosities. For each value of the integrated luminosity the signihcances are shown considering 
K = 10% systematic uncertainty. 


that the Hx distribution is extended well beyond 1 TeV. The optimized Hx cut of 800 GeV kills the ttHsM 
background by about 93%, while the tt -bjets background is entirely wiped out. After all cuts, the only 
surviving background turns out to be ttHsM as expected, with the background cross section being 0.028 fb. 

The cut efficiency after all the selection cuts, in the jet substructure scenario, increases steadily from the 
benchmark points PI to P5 due to a slight increase in the Higgs tagging efficiency as well as an increase in the 
strength of Hx and cuts. However the rapid fall of the stop pair production cross-section with increasing 
mass results in a minuscule final cross-section for stop masses above 1 TeV. To project the reach at 13 TeV, 
we define the signal significance (5) as S/\/(B -b (kB)^), where we include a systematic uncertainty factor 
of K in the background estimation to compensate for the fact that we do not perform any dedicated detector 
simulation in this study. For this study we choose k to be 10%. In Table 6, the signal significances for 100, 
300, 1000 fb~^ of integrated luminosity are tabulated in columns 6-8 for the five benchmark points. It can be 
noted that one can probe a light stop mass up to 800 GeV(Pl) at 100 fb“^ integrated luminosity at 13 TeV 
LHC with a significance > 5, while stop masses up to 1 TeV (P3) can be observed at 300 fb“^ luminosity. 
The signal suffers because of low production cross section, and only with the very high luminosity run at 
LHC (i.e 3000 fb“^), one can achieve signal significances greater than 5cr for the benchmark points P4 and 
P5 with stop masses greater than 1 TeV. We also observe that one could exclude mj ^ 1.2 TeV for 300 fb“^ 
integrated luminosity with the jet substructure technique. 

We compare the results obtained using the jet substructure method with the leptonic searches for this 
channel. We first analyze the case of single lepton as denoted by Case-1 in the previous section. In Table 
7 and Table 8, we summarize the event yields (normalized to cross section) for the signal and backgrounds 
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Effective 

cross-section after the cuts (in fb) 

Signal 

Production 

c.s. (fb) 

C6 

C7 

C8 

C9 

CIO 

Cll 

C12 

PI 

27.2 

16.39 

10.44 

3.55 

0.442 

0.441 

0.27 

0.239 

P2 

12 

8.49 

5.5 

1.88 

0.315 

0.314 

0.2 

0.176 

P3 

6 

4.34 

2.54 

0.88 

0.259 

0.259 

0.169 

0.139 

P4 

1.5 

1.17 

0.74 

0.246 

0.109 

0.109 

0.077 

0.061 

P5 

0.5 

0.45 

0.29 

0.098 

0.040 

0.040 

0.028 

0.0255 


Table 7: One lepton: Event summary for the signal after individual cuts as described in the text. In the 
last column we show the final cross-section (c.s.) after all the event selection cuts have been applied. 



Effective cross-section after the cuts (in fb) 

SM bkgs 

Production 

C6 

C7 

C8 

C9 

CIO 

Cll 

C12 


c.s. (fb) 








tf-l- jets 

700000 

20465.1 

8130.2 

1410.6 

9.13 

8.96 

5.6 

2.28 

ttHsM 

500 

54.75 

37.12 

8.1 

0.089 

0.089 

0.054 

0.036 

W -1- jets 

83315178 

6153.2 

0 

0 

0 

0 

0 

0 

Z + jets 

42070505 

4115.1 

470.3 

0 

0 

0 

0 

0 

Total 









bkg 








2.316 


Table 8: One lepton: Event summary for the backgrounds (bkgs) after individual cuts as described in the 
text. In the last column we show the final cross-section (c.s.) after all the event selection cuts have been 
applied. For the tt-f jets background we have generated a matched sample of tf-|- 0 jet, tt-l- 1 jet and tt-f 2 
jets using Madgraph. 



Signal(Ns) ( Background(NB)) 

Significance(<S) for 

K = 10% 


m5,(GeV) 

100 fb"i 

300 fb“i 

1000 fb”i 

100 fb"i 

300 fb"i 

1000 fb“i 

PI 

804.3 

23.9 (231.6) 

71.7 (694.8) 

239 (2316) 

0.86 

0.96 

1.01 

P2 

908.2 

17.6 (231.6) 

52.8 (694.8) 

176 (2316) 

0.64 

0.71 

0.74 

P3 

1003.7 

13.9 (231.6) 

41.7 (694.8) 

139 (2316) 

0.51 

0.56 

0.58 

P4 

1211.5 

6.1 (231.6) 

18.3 (694.8) 

61 (2316) 

0.22 

0.25 

0.26 

P5 

1392.5 

2.55 (231.6) 

7.65 (694.8) 

25.5 (2316) 

0.09 

0.10 

0.11 


Table 9: One lepton: The summary of our signal and backgrounds. Columns 3-5 show the number of signal 
(total background) events for three values of the integrated luminosity: 100 fb~^, 300 fb“^ and 1000 fb~^. 
The columns 6-8 show the statistical significance of our signal for the above three integrated luminosities. For 
each value of the integrated luminosity the significance is shown for three choices of the amount of possible 
systematic uncertainties, k = 10%. 


respectively. We find that the W-f jets background is entirely killed by the demand of b-tagged jets, as 
expected. For the Z -|- jets background, the demand of two b-tagged jets, along with four hard jets implies 
that no isolated lepton is expected to be present. Therefore, imposition of the single lepton criteria (see 
column C8, Table 8) eliminates the Z+ jets background. The remaining background consists of tt + jets 
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and ttHsM- As can be seen from column 9 in Table 8, plrp along with Ht and Mt suppresses this set to a 
negligible level. However the small starting cross section for the signal essentially leads to the fact that the 
signal signihcance is not large enough for a discovery (or an exclusion) even at the 13 TeV high luminosity 
scenario for these benchmark points as can be observed from Table 9. 

We next turn our attention to the di-lepton channel, denoted by case-II in the previous section. The 
signal and background yields in this case are tabulated in Table 10 and Table 11 respectively corresponding 
to the selection cuts C13 - C17 described in the last section. The demand for two isolated leptons in cut 
C13 suppresses the W -I- jets background, while rest of the background is suppressed by the Mt 2 cut (C16). 
The Mt 2 cut is expected to be bounded from above by the mass of the W boson for the background, while 
for signal it is expected to extend to higher values due to the presence of large missing energy. We observe 
that the backgrounds like tt -I- jets is reduced to zero (0) after the Mt 2 cut, while there is a small residue 
from the tfHsM background, which are further reduced by the cut of C17. In Table 12, we observe that 
for the di-lepton channel, up to 1 TeV stops can be discovered with 100 fb~^ integrated luminosity, while a 
light stop mass up to 1.2 TeV can be discovered with 300 fb“^ of luminosity. Therefore we conclude that 
the di-lepton channel is most suited to probe the heavier ti mass with the benchmark points depicted in 
Table 2. We also observe that stop masses up to 1.4 TeV can be excluded at 13 TeV LHC, with 300 fb“^ 
luminosity. 




Effective cross- 

section after the cuts (in fb) 

Signal 

Production 

c.s. (fb) 

C13 

C14 

C15 

C16 

C17 

PI 

27.2 

24.61 

15.29 

1.49 

0.17 

0.038 

P2 

12 

11.31 

7.17 

0.63 

0.096 

0.0264 

P3 

6 

5.65 

3.19 

0.297 

0.061 

0.0262 

P4 

1.5 

1.43 

0.875 

0.065 

0.019 

0.012 

P5 

0.5 

0.494 

0.314 

0.026 

0.0069 

0.0035 


Table 10; Two lepton: Event summary for the signal after individual cuts as described in the text. In the 
last column we show the final cross-section (c.s.) after all the event selection cuts have been applied. 



Effective 

cross-section after the cuts (in fb) 

SM bkgs 

Production 

C13 

C14 

C15 

C16 

C17 


c.s. (fb) 






tf-t- jets 

700000 

105014.18 

38630.79 

515.89 

0 

0 

ttHsM 

500 

179.56 

116.06 

3.01 

0.007 

0.001 

W -1- jets 

83315178 

123064.92 

0 

0 

0 

0 

Z -1-jets 

42070505 

78891.62 

4799 

0 

0 

0 

Total 







bkg 






0.001 


Table 11: Two lepton: Event summary for the backgrounds (bkgs) after individual cuts as described in 
the text. In the last column we show the final cross-section (c.s.) after all the event selection cuts have been 
applied. For the tt-l- jets background we have generated a matched sample of tt-l- 0 jet, tt-l- 1 jet and tt-l- 2 
jets using Madgraph. 


Finally, we examine the possibility of the detection of multiple Higgs peaks in the fat-jet mass (mj) 
distribution over the SM background, where J denotes the fat-jet. In Fig. 5, we plot the fat-jet distribution 
for the jet substructure technique at 300, 1000, and 3000 fb“^ luminosity respectively. We plot first six 
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Signal(Ns) ( Background(NB)) 

Significance (5) for 

K = 10% 


mq (GeV) 

100 fb"i 

300 fb“i 

1000 fb“i 

100 fb"i 

300 fb"i 

1000 fb"i 

PI 

804.3 

3.8 (0.1) 

11.4 (0.3) 

38 (1) 

12.1 

20.8 

37.8 

P2 

908.2 

2.64 (0.1) 

7.92 (0.3) 

26.4 (1) 

8.34 

14.4 

26.3 

P3 

1003.7 

2.62 (0.1) 

7.86 (0.3) 

26.2 (1) 

8.28 

14.3 

26.0 

P4 

1211.5 

1.2 (0.1) 

3.6 (0.3) 

12.0 (1) 

3.79 

6.56 

11.9 

P5 

1392.5 

0.35 (0.1) 

1.1 (0.3) 

3.5 (1) 

1.1 

2.0 

3.48 


Table 12: Two lepton: The summary of our signal and backgrounds. Columns 3-5 show the number 
of signal (total background) events for three values of the integrated luminosity: 100 fb“^, 300 fb“^ and 
1000 fb“^. The columns 6-8 show the statistical significance of our signal for the above three integrated 
luminosities. For each value of the integrated luminosity the significance is shown for three choices of the 
amount of possible systematic uncertainties, k = 10%. 


PT-ordered fat-jet masses at the end of our selection cut (C5) in our jet substructure Higgs analysis. In 
these plots we demand two b-tagged jets only for the fat-jet satisfying the 125 GeV Higgs criteria, i.e, 
110 < mj < 140 GeV. We observe that while the SM Higgs peak is overwhelmingly visible at all the 
luminosity options, one can also observe slight excesses in three different masses over the ttHsM background. 
We find small excesses in the lighter CP even Higgs boson (~ 65 GeV), Z boson (~ 91 GeV) as well as at the 
top quark (~ 172 GeV) masses. With a better branching ratio of x)* —t Xi^i {i = 2,.., 5), it is possible that 
the peak corresponding to the lighter CP even Higgs might be visible at even lower luminosities. However, 
for the benchmark points of our choice which are not the most ideal for this purpose as the branching ratio 
X 2 XiHi is not large enough, this possibility can only be realized at very high luminosity options. Using 
a simple number counting method we observe that with 3000 fb“^ luminosity one can achieve a 3cr signal 
signihcance to observe the light CP-even Higgs boson {Hi) for our representative benchmark point P3 with 
stop mass around 1 TeV. The possibility of having significant BR for the decay of X 2 3 to both the Higgses 
(Hi, sm) is beyond the scope of this paper, and we try to report this possibility elsewhere [73]. 



nij (GeV) 


Figure 5: Fat-jet mass distribution with 300 (left), 1000 (middle), 3000 (right) fb“^ luminosity, after all 
cuts as described by the jet substructure analysis in the text. Here we plot first six px-ordered fat-jet masses 
at the end of our selection cut (C5), no extra b-tagging criteria on these fat-jets are imposed expect the SM 
Higgs like fat-jets where we demand two b-tagged jets inside the fat-jet. The reason behind this plot is to 
investigate the possibility of having multiple Higgs peaks in the fat-jet mass distributions. The figures are 
all normalized to cross section with a particular choice of luminosity. 
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5 Summary 

In this article, we studied Higgs signatures from the cascade decay of the light stop (ti) in the framework of 
NMSSM. We identify regions of parameter space where this channel is dominant and consistent with current 
Higgs coupling measurements, LHC constraints, dark matter relic density and direct detection constraints, 
and provide a detailed collider simulation of the signal and background. We probe this channel using the jet 
substructure techniques and via the conventional leptonic searches. It is concluded that in terms of signal to 
background ratio, the di-lepton channel is the most promising one. In the di-lepton channel we can discover 
stop masses up to 1.2 TeV with 300 fb“^ luminosity. With the jet substructure method, stop masses up to 1 
TeV could be probed with 300 fb“^ luminosity. Finally we also investigate the possibility of the appearance 
of multiple Higgs peaks in the fat-jet mass distribution. This is plagued by the low branching ratio to the 
lighter CP even state in our case, and is only visible at very high luminosity run of the LHC. 
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